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ABSTRACT 


Numerical  values  of  static  stress  concentration  factor  for  ten 
typical  tunnel  cross-sections,  determined  by  the  methods  of  photo¬ 
elasticity  are  presented  herein.  Results  indicate  that  stress 
distribution  around  a  cavity  are  very  sensitive  to  its  cross-section’s 
configuration  and  the  direction  of  the  applied  load.  The  values  of 
the  stress  concentration  factors  for  the  configuration  considered 
ranges  from  3.24  -  6.10  for  maximum  compressive  stress  and  ranged 
from  1.07  -  2.18  for  the  maximum  tensile  stress. 
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1.  OBJECTIVE 

It  was  desired  to  obtain  the  elastic  solution  of  boundary  stresses 
developed  in  deep  tunnels  under  uniform  uniaxial  load.  The  study  is 
limited  to  elastic,  homogeneous  materials.  Ten  prescribed  tunnel  con¬ 
figurations  were  investigated. 

2 .  INTRODUCTION 

Recent  studies  concerning  the  feasibility  of  deep  underground 
command  and  control  installations  have  indicated  that  a  cavity  located 
under  three  to  four  thousand  feet  of  overburden  could  withstand  the 
shock  due  to  the  direct  hit  of  a  multi-megaton  weapon.  This  conclusion 
is  based  on  the  extrapolation  to  depth  of  the  free-field  stresses  in 
the  close-in  region  where  a  hydrodynamic  analysis  is  valid. However, 
the  survivability  of  a  deep  underground  cavity  depends  not  only  on  the 
magnitude  of  the  free-field  stress  but  also  on  geometric  configuration 

of  the  cavity  itself. 

(2) 

Baron,  et  al  have  investigated  the  diffraction  effects  of  a 
circular  cavity  in  an  elastic  half-space,  and  have  shown  that  there  is 
a  maximum  dynamic  stress  magnification  of  approximately  3.3;  the 

(3) 

corresponding  static  stress  magnification  being  3.,  ARF  has  conducted 
a  photo-elastic  study  where  the  dynamic  stress  concentration  factors 
for  circular,  elliptical  and  square  (with  filleted  corners)  holes  are 
determined.  It  was  found  that  the  difference  between  the  maximum  static 
and  maximum  dynamic  stress  is  of  the  order  of  10-20%. 

In  view  of  this  apparent  static -dynamic  similitude  and  since  the 
static  photo -elastic  analyses  are  much  easier  to  carry  out,  it  was  felt 
that  a  static  photo-elastic  investigation  of  the  boundary  stresses 
associated  with  several  typical  cavity  cross-sections  was  warranted. 

This  report  presents  the  results  of  such  an  investigation. 

3.  EXPERIMENTAL  APPROACH 

The  elasticity  problem  of  a  long  tunnel  of  constant  cross-section, 
subjected  to  uniform  loading  along  its  length  is  in  the  category  of 
plane  strain.  It  is  shown  in  the  Theory  of  Elasticity  that  the  stress 
distribution  developed  in  planes  of  the  cross-section,  i.e.,  in  planes 
normal  to  the  axis  of  the  tunnel,  is  identical  to  that  developed  in 
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the  equivalent  plane-stress  system.  Thus,  thin  models  were  employed 
in  which  the  tunnel  dimensions  were  large  compared  to  the  model  thick¬ 
ness. 

While  the  loading  system  of  interest  is  uniaxial  compression,  the 
system  employed  in  these  tests  was  uniaxial  tension.  It  is  clear  that 
the  distribution  of  stresses  in  the  small  strain  elasticity  solution  is 
independent  of  the  sign  of  applied  load.  The  boundary  stresses  shown' 
in  the  accompanying  graphs  are  for  tensile  loads;  to  convert  to  the 
equivalent  system  of  compressive  loads,  interchange  symbols  for  tensile 
and  compression  stress. 

Ideally,  a  deep  tunnel  in  a  continuous  medium  should  be  represented 
by  a  model  whose  overall  width  is  large  compared  to  the  tunnel  width. 

As  a  practical  compromise,  a  ratio  of  4:1  was  chosen;  with  these  pro¬ 
portions,  the  proximity  of  the  free  edges  of  the  model  have  little 
influence  on  the  stress  distribution  on  the  tunnel  boundary. 

The  photoelastic  method  was  chosen  for  the  investigation  of  these 
irregular  tunnel  cross-sectional  shapes,  rather  than  the  analytical 
method  of  elasticity.  Since  the  boundary  shapes  could  be  produced 
mechanically  but  could  not  be  satisfactorily  prescribed  by  a  mathematical 
equation,  the  experimental  method  was  by  far  the  easier  of  the  two 
alternatives . 

Photoelasticity  is  based  upon  a  unique  property  of  certain  trans¬ 
parent  materials,  notably  some  plastics.  When  a  ray  of  light  enters  a 
stressed  photoelastic  model,  it  divides  into  two  plane  polarized  components, 
with  directions  of  polarization  parallel  to  the  directions  of  principal 
stresses.  In  addition,  each  component  propagates  through  the  model 
material  with  a  velocity  that  depends  upon  the  state  of  stress,  and  the 
two  components  emerge  out  of  phase  with  each  other.  The  components  are 
recombined  in  the  photoelastic  polariscope  to  exhibit  constructive  or 
destructive  interference,  according  to  the  phase  difference.  In  this 
way,  a  continuous  pattern  of  optical  interference  fringes  (or  bands)  is 
produced  throughout  the  field  of  view.  The  phase  difference  at  any 
point  in  the  model  can  be  determined  from  the  interference  pattern,  and 
from  this,  information  relating  to  the  state  of  stress  at  any  point  is 
obtained. 
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The  phase  difference,  or  in  more  common  terminology,  the  photo¬ 
elastic  fringe  order,  is  proportional  to  the  difference  of  principal 
stresses  developed  at  each  point  in  the  model.  At  free  boundary  points, 
e.g.,  the  tunnel  boundaries  in  these  analyses,  the  principal  stress  normal 
to  the  boundary  must  be  zero,  and  the  fringe  order  is  directly  proportional 
to  the  tangential  boundary  stress.  Thus,  boundary  stress  distribution  is 
given  by  the  graph  of  fringe  order  vs.  position  along  the  boundary. 

4.  MODEL  AND  LOADING  SYSTEM 

A  commercial  plastic,  Homolite  100,  was  selected  as  the  model 
material.  It  exhibits  high  stress-optical  sensitivity  and  other  desirable 
photoelastic  properties.  It  is  available  in  large  sheets  as  a  colorless, 
homogeneous  plastic  with  glass-like  surfaces. 

Model  dimensions  were  10"  width,  40n  length,  3/16"  thickness,  with 
2  1/2M  tunnel  width.  Tunnel  configurations  are  given  in  Table  I.  Loads 
were  applied  by  means  of  four-arm  whipple-trees  as  illustrated  in  Fig.  1. 

Well  machined  whipple-trees  were  employed  to  produce  a  uniform  load  across 
the  plate  without  introducing  any  transverse  constraints;  their  effectiveness 
is  evidenced  by  the  striking  symmetry  of  the  photoelastic  patterns  shown  in 
succeeding  figures. 

The  applied  tensile  loads  were  not  measured.  Instead,  a  self- indicating 
technique  was  used  in  which  measurements  were  taken  in  a  region  where  the 
relationship  between  stress  and  load  was  the  same  for  all  specimens.  For 
increased  accuracy,  a  stress  concentration  in  the  form  of  a  small  hole 
(3/4n  diameter)  was  introduced  midway  between  the  tunnel  and  whipple- 
tree  and  the  maximum  fringe  order  at  the  hole  boundary  was  observed  in 
each  test.  This  data  permitted  direct  correlation  between  models  on  the 
basis  of  a  common  tensile  load. 

5.  PHOTOELASTIC  POLARISCOPE 

A  diffuse  light  doubling  polariscope  was  employed,  as  illustrated  in 
Fig,  1.  In  this  system,  a  mirror  placed  immediately  behind  the  model  returns 
the  light  rays  for  a  second  passage  through  the  model,  thus  doubling  the  phase 
difference  and  doubling  the  sensitivity  of  the  experiment.  Light-field 
circular  polarization  was  employed;  in  order  to  eliminate  the  influence  of 
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polarization  by  the  beam-splitter,  the  polarizer  axis  was  vertical,  the 
analyzer  axis  horizontal,  and  the  quarter -wave  plate  axes  45°  from  the 
vertical.  A  sodium  vapor  lamp  was  used  as  the  monochromatic  light  source 
The  camera  was  equipped  with  a  20"  focal  length  lens  and  a  Polaroid-Land 
camera  back;  Polaroid  film  yielding  positive  transparencies  was  used. 

The  entire  polariscope  was  mounted  on  an  elevator  pedestal  so  that 
photographs  of  the  tunnel  and  reference  hole  could  be  taken  in  sequence. 

6.  EXPERIMENTAL  PROCEDURE 

The  optical  system  was  aligned  by  optical  reflection  techniques  to 
position  the  planes  of  the  model  and  mirror  perpendicular  to  the  axis  of 
the  camera.  The  camera  was  focused  on  the  plane  of  the  model. 

For  each  tunnel  shape,  the  photoelastic  model  was  machined  by  means 
of  a  high  speed  routing  machine  equipped  to  follow  a  previously  prepared 
metal  template  of  the  model  contours.  Tungsten  carbide  tools  were  used. 

A  jet  of  compress  air  was  played  on  the  tool  for  cooling  and  chip  removal 

The  model  was  cleaned  and  installed  in  the  loading  frame  without 
delay.  A  deleterious  time-edge  effect  gradually  develops  in  the  model 
material  upon  exposure  of  freshly  cut  surfaces  to  the  atmoshere,  but  this 
remains  inconsequential  if  the  experiment  is  performed  rapidly.  The 
tensile  load  was  increased  until  the  photoelastic  fringe  order  at  the 
most  highly  stressed  point  was  approximately  twelve. 

The  photoelastic  pattern  surrounding  the  tunnel  was  photographed. 

The  optical  system  was  then  raised  to  the  level  of  the  reference  hole 
and  this  photoelastic  pattern  was  photographed. 

Positive  photographic  prints  were  produced  for  each  tunnel  pattern 
as  shown  in  Figs.  2  through  11. 

The  fringe  orders  along  the  tunnel  boundaries  were  ascertained, 

and  divided  in  each  case  by  the  maximum  fringe  order  at  the  reference 

hole.  This  dimensionless  ratio  was  plotted  to  such  a  scale  that  the 

maximum  stress  developed  in  the  case  of  the  circular  tunnel,  Fig.  2, 

would  equal  3,24a  .  Here,  a  represents  the  tensile  load  divided  by 
o  o 

the  gross  cross-sectional  area  (10"  x  3/16")  of  the  plate;  alternatively, 
<jQ  represents  the  tensile  stress  that  would  have  existed  in  the  region 
of  the  tunnel  if  the  tunnel  was  not  present.  The  factor  3.24  is  the 
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known  stress  concentration  factor  for  a  circular  hole  in  a  finite 
tensile  plate  of  these  proportions,  when  based  upon  the  gross  cross- 
sectional  area.  Thus,  the  graphs  shown  in  Figs.  11  through  21  are 
unified  on  the  basis  of  a  common  applied  stress. 

7  .  RESULTS 

The  results  of  the  investigation  are  presented  in  graphical  form 
in  Figs.  12  through  21.  In  each  case  the  magnitude  of  stress  at  any 
boundary  point  is  given  by  the  length  of  the  generatrix  normal  to  the 
boundary . 

Peak  tensile  and  compressive  stresses  are  given  in  numerical  form 
on  the  graphs  and  also  in  Table  II. 

The  limit  of  experimental  error  is  estimated  as  *1/4  of  a  fringe 
order.  Maximum  fringe  order  at  the  reference  hole  averaged  7  for  all 
tests,  yielding  a  representative  maximum  error  of  approximately  *4%  in 
unifying  the  data  on  the  basis  of  a  common  load.  Maximum  fringe  order 
at  the  tunnel  boundary  averaged  approximately  10,  yielding  a  representative 
maximum  error  of  *2  1 /2%  of  the  maximum  stress.  Accordingly,  the  accuracy 
of  individual  stress  distributions  shown  in  the  graphs  is  estimated  to  be 
within  1"2  l/27o  of  peak  stress,  while  the  accuracy  of  corresponding  numerical 
values  is  within  *6  1/2%. 

The  results  represent  the  effect  of  uniaxial  loading  applied  in  the 
direction  of  the  vertical  axis  of  symmetry  for  tunnel  Models  1  through  9. 

In  the  case  of  Model  10,  the  uniaxial  loading  was  applied  at  an  angle  of 
30°  from  the  axis  of  symmetry. 

8.  CONCLUSION 

Observing  the  result  obtained  herein,  one  would  recommend  strongly 
that,  tunnel  shapes  shown  in  Figures  12,  13,  13,  16,  18,  and  20  be  used 
for  any  future  underground  command  and  control  installation,  or  in  general 
it  is  preferred  that  the  axis  of  the  tunnel  which  lies  in  the  direction  of 
loading  be  larger  than  the  axis  that  is  perpendicular  to  the  line  of  loading. 

Another  very  interesting  and  perhaps  very  important  result  obtained  in 
this  report  is  the  extreme  sensitivity  of  the  stresses  around  the  cavity 
due  to  a  change  of  direction  of  loading.  This  is  shown  in  Figure  21  which 

a 

Peterson,  R.  E.,  Stress  Concentration  Design  Factors,  John  Wiley  and  Sons,  Inc 6 , 

New  York,  1953. 


has  the  same  cross-section  as  Figure  13  except  the  load  is  applied  at  30° 
to  its  axis  of  symmetry.  In  the  case  of  21,  it  is  found  that  maximum 
compressive  stress  is  5.10  and  maximum  tensile  stress  is  2.18  as  compared 
with  3.24  and  1.26  in  Figure  13.  Therefore,  great  care  must  be  exercised 
in  choosing  the  tunnel  configuration  so  that  all  factors  are  to  be 
incorporated.  Ideally,  one  should  choose  the  circular  cross-section  so 
that  the  axis  of  symmetry  will  be  preserved  in  all  directions. 
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FIG. 18 
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FIG.  21 

MODEL  10 


TABLE  I 


TUNNE L  CONFIGURATIONS 


Model  No 

1 

2 

3 

4 
3 
6 

7 

8 
9 

10 


h  height  of  rectangle 
r  radius  of  circle 
r^  major  radius  of  ellipse 
r^  minor  radius  of  ellipse 


Description 

Circle 
3/4  Circle 
Semi-circle 

Semi-circle  above  rectangle,  h  *  r 
Semi-circle  above  square 
Shallow  semi-ellipse,  r^/^  “  1*5 
Shallow  ellipse  above  rectangle,  h  ■  r^ 

Deep  semi-ellipse,  r^/r ^  =  1*5 
3/4  Ellipse  (cut  at  1.5r^) 

3/4  circle  oriented  with  axis  of  symmetry 
30°  from  direction  of  load 


Note:  Tunnel  width  held  constant  at  2  1/2”;  1/8"  radius  fillets 

at  all  corners. 
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TABLE  II 


maximum  stresses 


Spec .  No . 

Max.  Tensile 

Max.  Compressive 

Stress 

Stress 

X 

3.24a 

o 

1 . 36a 

o 

2 

3.24a 

o 

1.26a 

o 

3 

5.95a 

o 

1.19a 

o 

4 

4.07a 

0 

1 . 20a 

o 

5 

3.90a 

o 

1.43a 

o 

6 

6.10a 

o 

1.25a 

0 

7 

4.20a 

o 

1.07a 

o 

8 

4.80a 

o 

1.20a 

o 

9 

2.69a 

o 

1.34a 

o 

10 

3.10a 

o 

2.18a 

o 
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